Entanglement negativity in quantum field theory 
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We develop a systematic method to extract the negativity in the ground state of a 1+1 dimensional 
relativistic quantum field theory, using a path integral formalism to construct the partial transpose 
p'^ of the reduced density matrix of a subsystem A — Ai U A2, and introducing a replica approach 
to obtain its trace norm which gives the logarithmic negativity £ = ln\\p'^\\. This is shown to 
reproduce standard results for a pure state. We then apply this method to conformal field theories, 
deriving the resuh £ ~ (c/4) In (^1^2/(^1 +^2)) for the case of two adjacent intervals of lengths £i,£2 
in an infinite system, where c is the central charge. For two disjoint intervals it depends only on the 
harmonic ratio of the four end points and so is manifestly scale invariant. We check our findings 
against exact numerical results in the harmonic chain. 
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Recent years have witnessed a large effort to under- 
stand and quantify the entanglement content of many- 
body quantum systems (see [T] for reviews). This is 
usually achieved by partitioning an extended quantum 
system into two complementary subsystems and calcu- 
lating the entanglement entropy Sa, defined as the von 
Neumann entropy of the reduced density matrix pA of 
one subsystem. However, this procedure does not give 
information about the entanglement between two non- 
complementary parts Ai and A2 of a larger system be- 
cause generically their union is in a mixed state. The 
mutual information Sai + Sa2 ^ SaiuA2 measures the 
correlations between the two parts, but gives only an up- 
per bound on the entanglement between them. 

A more useful measure of entanglement in this case is 

the negativity pi, defined as follows. Denoting by \e[^^) 
(2) 

and \ej ) two bases in the Hilbert spaces Hi and H2 of 
each part, one first defines the partial transpose of p as 
(ef^efV^^lei'^ep)) = {e^'>4'^\p\e^^^ef) and then the 
logarithmic negativity as 

£ = ln||p^^|| = InTrlp'^^^l , (1) 

where the trace norm is the sum of the absolute 

values of the eigenvalues of p^^ . When the two parts 
are two microscopic degrees of freedom (e.g. spins), the 
negativity coincides with other commonly used entangle- 
ment estimators [U |3] , but its definition is more appeal- 
ing because it is basis independent and so calculable by 
quantum field theory (QFT). 

The use of QFT naturally unveils universal features, in 
particular close to a quantum critical point. For ID crit- 
ical theories, that at low energy are also Lorentz invari- 
ant, the powerful tools of conformal field theory (CFT) 
can be applied. As a matter of fact, the interest in en- 
tanglement in extended systems has been considerably 
boosted by the now classical CFT result that the entan- 



glement entropy of a large block of length £ is Sa — ^ In £, 
with c the central charge [IHS]- When a subsystem con- 
sists of two blocks, the entanglement entropy can also 
be obtained from CFT [31^1, but this gives only the mu- 
tual information, not the entanglement between the two 
blocks. 

For these reasons, and also motivated by recent results 
in some ID models [TDHH], in this Letter we carry out 
a systematic study of the logarithmic negativity in QFT 
(in particular CFT) based on a new replica formalism. 

A replica approach. We consider the traces of integer 
powers of p"^^. For n even (odd), let us say rie — 2m 
{rio = 2m +1), they read 

Tr(p^^)"= = ^ A^^ = ^ \kr^ + ^ Ihr- , (2) 

i li>0 li<0 

Tr(p^^)"° = E^"° = Ei^^r-Ei'«r- 

i li>0 li<0 

The analytic continuations from even and odd n are dif- 
ferent and the trace norm in which we are interested is 
obtained by considering the analytic continuation of the 
even sequence at rig — >■ 1, i.e. £ = lim lnTr(p^^)"'=, 

while the limit rio — !■ 1 gives the normalization Trp-^^ = 1. 

As a first example, let us consider the case in which 
p — corresponds to a pure state \'^). Then, the 

eigenvalues of p^^ are related to the Schmidt decomposi- 
tion coefficients [H [T3] and after simple algebra 

Tr(/=^)"= = (Trp^=/^)2 , Tr(p^^)"° = Trp^ , (3) 
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FIG. 1: We consider the entanglement between two blocks Ai 
and A2 embedded in the ground-state of a larger system. 



FIG. 2: Top: The reduced density matrix pA of two disjoint 
intervals. Middle: Partial transpose with respect to the sec- 
ond interval . Bottom: Reversed partial transpose . 

where p2 is the reduced density matrix on 'H2- Taking 
the limit rig — >■ 1, we recover the result 2J that for a 
pure state the logarithmic negativity is the Renyi entropy 
Si/2 = 21nTrp2^^ 

Negativity and QFT. For concreteness we refer to a ID 
system and we consider the tripartition depicted in Fig. [l] 
with A composed of two parts A — AiU A2 ^ [ui,vi] U 
[u2,V2] a-nd B the remainder, but most of the following 
ideas apply to more general cases. In the ground-state 
of a QFT, the reduced density matrix pA has the path 
integral representation in Fig. [2] (top) [5] . The two open 
cuts correspond to the rows and columns of pA- Trp^ for 
integer n can be obtained by joining cyclically n of the 
above density matrices as in Fig. |3] (top). Thus Trp^ is 
(proportional to) the partition function on this n-sheeted 
Riemann surface which is equivalent to the correlation 
function of the twist fields Tn{z) constructed exploiting 
the cyclic permutation symmetry of the sheets, i.e. [6l [7] 

Trp:^ = {Tn{u,)fnMTn{u2)fn{v2)) ■ (4) 

The partial transposition with respect to the second in- 
terval A2 corresponds to the exchange of row and column 
indices in A2. In the path integral representation, this 
is equivalent to interchange the upper and lower edges 
of the second cut in pA as in the middle of Fig. [2] It is 
convenient to reverse the order of the column and row 

indices in A2 as in the bottom of Fig. [21 to obtain the 

c • • — 

reversed partial transpose p^^. This is related to the par- 
tial transpose as p^^ — Cp^C, where C reverses the 
order of indices either on the lower or on the upper cut. 
Clearly Tr(p^^)" = Tr(p5')" and so Tr(p^^)" is the parti- 
tion function on the n-sheeted surface obtained by joining 
cyclically n of the above p^^ as in the bottom of Fig. ^ 
It is then straightforward to see that 

T^^PaT = {Tn{ui)Tn{vi)fn(u2)Tn{v2)) , (5) 

i.e. the partial transposition has the net effect to ex- 
change two twist operators compared to Eq. ^^. To 
replace p^ with pj^ it has been fundamental to consider 
integer cyclical traces. The operator C enters in quanti- 
ties like Tr(p^p^^) which is in fact the partition function 
on a non-orientable surface with the topology of a Klein 
bottle. This can be computed using CFT methods [14] . 




FIG. 3: Path integral representation of TrpJ (top) and 
Tr(p^^)" (bottom) for n = 3. 

For n ^ 2, T2 ^ T2 and so Trp^ = Tr(pJ)2 which 
follows from the properties of the trace. 

We first specialize to a pure state by letting S for 
which Tr(p^^)" can be worked out in full generality as 

Triple ^{T^iu2)T:iiv2)). (6) 

This expression depends on the parity of n because 
connects the j-th sheet with the {j + 2)-th one. For n = 
rig even, the rie-sheeted Riemann surface decouples in 
two independent (ne/2)-sheeted surfaces. Conversely for 
n = rio odd, the surface remains a rio-sheeted Riemann 
surface. Thus we have 

Tr(pJ)"= = {{%.j2{u2)%j2{v2))f = WA!'f^ 

Tr(p5^)"° - {TnM2)%Av2)) = Trp:^; , (7) 

which are the results for pure states in Eq. ([s]), recovered 
here purely from QFT. 

We now specialize to the case of a CFT, for which the 
twist fields transform like primary operators of dimension 
At"^ — c{n — l/n)/12 [B]. Thus when A2 is embedded in 
an infinite system we have {i — U2 — V2) 

Tr(pJ)"= ocrt^'^-^^), Tr(pJ)"° cxrs'""-^). (8) 

Despite of the simplicity of the above calculation, it 
shows one important point of the CFT analysis: for 
n = rif, even, T^^^ has dimension A7-2 = c(ne/2— 2/ne)/6, 
while for n = rio odd, T^^ has dimension A7-2 = 
c{no — l/7io)/12, the same as Tn^- We finally have 

IIp^^II = Jirn^Tr(p^^)"= cx£5^ f^^ln^-f cnst . (9) 
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Two adjacent intervals. Let us now consider the non- 
trivial configuration in which two intervals Ai and A2 of 
length £1 and ^2 share a common boundary (let us say at 
the origin) which is described by the 3-point function 



(10) 



whose form is determined by conformal symmetry |15j . 
For n — Uf. even, using the dimensions of the twist oper- 
ators calculated above, we find 

Tr(pJ)"= oc (44)"'^^"^'(4 +^2)"^^'^ + ^^ , (11) 
that in the limit — > 1 gives 



T2\ 



^5^^1n^+cnst. (12) 



For n^Uo odd, Tr(pJ)"° cx (^i4(^i -f 4)) 



:) 



that for 



1 gives again Trpy 



1. 



All the previous results may be generalized to the case 
of a finite system by using a conformal mapping from the 
cylinder to the plane. This results in replacing I with the 
chord length {L/n) sin(7r^/L). 

Two disjoint intervals. For the more interesting and 
complicated situation of two disjoint intervals of Fig. [l] 
global conformal invariance gives {li = \vi — Ui\) 

Tri/^r « [4^2(1 - 2;)]-^("-^^t/n(2/) , (13) 

where y = (•"i~"i)(^2-"2) ■ ^j^g four-point ratio (0 < u < 
1) and Qn{y) a function depending on the full operator 
content of the theory. Trp^ in Eq. Q admits the same 
scaling form, but with a different scaling function J^n{y) 
which has been calculated for the free compactified boson 
and for the Ising model [THS]. Since Eqs. Q and ([s]) 
are related by an exchange of two twist fields, these two 
functions are related as 



^„(y)^(l-y)S("-^)j-„(2;/(y-l)) 



(14) 



Taking the replica limit rig — 1, we obtain 
S{y) = hm lng„^(y) = lim In [T^^ {y/{y - 1))] . (15) 

Then for conformal invariant systems, the negativity is a 
scale invariant quantity (i.e. a function only of y) because 
all the dimensional pref actors cancel in the replica limit. 
This has been argued already in the literature on the 
basis of numerical data [inilll], but never proved. 

In Refs. [3 [5] the function J-"„(x) has been obtained 
for some CFTs only for < x < 1 and it is a non-trivial 
technical problem to extend it to the domain a; < in 
which we are now interested. It is a hard open problem 
to find the analytic continuation to rig 1. We will 
report these technicalities for few specific cases elsewhere 
|16j and we limit here to discussing the main physical 



consequences of Eqs. (13), (14), and (15). These are 
highlighted by considering the limit y — )■ 1 and y — 0, i.e. 
close and far intervals respectively. If 1*2 ^ vi then y — > 1 
and we should recover the previous result for adjacent 



intervals. Comparing Eqs. ( 11 ) and ( 13 ) we have Gn{y) oc 
(1 — y)" (apart from possible multiplicative logarithmic 
corrections) with a equal to Aj-2 the dimension of 7^, i.e. 
ckrie = c(ne/2 — 2/ne)/6 and = c(no — l/no)/12. For 
TT-e — > 1 we have an^^i = i.e. the scaling function 

diverges approaching y = 1. The opposite limit of far 
intervals y — is worked out from the small y expansion 
of J-n (y) carried out in full generality in Ref. [5] . This is a 
sum over all intermediate operators of the form J-n{y) — 
J2iy'^^'^n{i)- The coefficients s„(i) have been explicitly 
calculated [5] and they do not depend on the parity of 
n. Thus, in the limit n — > 1 all these coefficients vanish, 
because the analytic continuation for even and odd n is 
the same (as the direct computation shows) and £{y) 
vanishes in y = faster than any power. 

The harmonic chain. We check the CFT results 
against exact computations in the harmonic chain with 
Hamiltonian 



1 



L 



(16) 



and periodic boundary conditions. For w = the chain 
is critical and its continuum limit is the c = 1 free boson. 
The construction of the partial transpose is detailed in 
[T7] and here we limit to presenting numerical checks of 
our CFT predictions. For w = 0, the zero mode leads to 
divergent expressions, thus we work at finite but small to 
such that ujL ^ 1. 

We first consider the case of two adjacent intervals of 
equal length £. The results for Tr(p^^)" for n = 3,4, as 
well as the results for the negativity £ are reported in 
Fig. |4] where they are compared with the finite size CFT 
predictions finding excellent agreement. 

The negativity of two disjoint intervals has been al- 
ready considered numerically |llj . We consider here the 
ratio 



(17) 



in which the non-universal parts due to the zero mode 
cancel and we are left with a universal function of y. 
The CFT prediction for this ratio is [TB] 



jl- y)l(-i'>UlZlFKiy)F^{l-y) 

nri!Re(j^,(^)F,(^)) 



(18) 

where Fq(x) =2Fi[q,l — q,\,x) being 2Fi{a,b,c, z) the 
hypergeometric function. This prediction is compared 
to the numerical data in Fig. [5j As i increases, the 
data approach the CFT result. The differences with the 



4 







■ L=150 


a)L=10"* 


• L=100 


00 L=10 " 


' L=150 


(0 L=10 ~ 




0.9- 



0.8 




0.02 0.04 0.06 0.08 

-2/3 





n=3 












\ % 










+ 






FIG. 4: For two adjacent intervals of equal length £ < L/2, 

we plot r„ = ln[Tr(p^'^^=')"/Tr(/9^-*=='""')"] as function of 
z = £/L. The subtraction is chosen to cancel non-universal 
factors. The bottommost panel shows £i = £ — (lnL)/4 in 
which non-universal terms are absent. The continuous lines 
are the parameter free CFT predictions. 



FIG. 5: Top: The ratio R„{y) in Eq. (17 1 as function of 



y for several L and for n = 3,4. The continuous lines are 
the parameter free CFT predictions. The inset shows a finite 
size scaling analysis for d„ = R'^^^{y) — Rnijj) for n = 3 
displaying the unusual correction [18]. The same is true 

for higher n |I6|. Bottom: The negativity £{y) is a universal 
scale invariant function with an essential singularity at t/ = 0. 



asymptotic formula are due to the presence of unusual 
corrections to the scaling [T5] of the form L~^/". A quan- 
titative finite size scaling analysis is reported in the inset 
of the figure. The bottom panel of Fig. [5] shows the neg- 
ativity £ for which all data collapse on a single curve, 
without sizable corrections. For small y, the data are 
very close to zero and are consistent with the form e""/** 
[llj . vanishing faster than any power. For y — > 1, we find 
e^^v) - (1 - y)-^'^\ ln(l - y)\''^/^ as obtained from the 
analytic continuation of Eq. ( [Tsl in this regime [TB] . The 
logarithmic correction may be responsible for the expo- 
nent ^ found in Ref. [llj as compared with our analytic 
result i, which is consistent with our general result |. 

Conclusions. We described a general QFT formalism 
to calculate the logarithmic negativity. For a conformal 
invariant theory we worked it out for two intervals, both 
adjacent and disjoint. In the latter case, the negativity is 
a universal scale invariant function. Some generalizations 
such as for compactified free boson, Ising CFT, finite 
temperature CFT, and massive QFT have been already 
obtained and will be presented elsewhere [IB] . 

However, there are still open problems, among them 
the analytic continuation rig —> 1 of the results for dis- 
joint intervals which remains a formidable task, reflecting 
a similar problem for the entanglement entropy [71 [H| . 

Finally, it is of extreme interest to check numerically 
our CFT predictions in more complicated lattice models 
such as spin-chains and itinerant fermions. 
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